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Stable Multiscale Petrov-Galerkin Finite Element 
Method for High Frequency Acoustic Scattering 

D. Gallistl* D. Peterseim* 


Abstract 

We present and analyze a pollution-free Petrov-Galerkin multiscale 
finite element method for the Helmholtz problem with large wave num¬ 
ber K as a variant of [Peterseim, arXiv:1411.1944, 2014]. We use stan¬ 
dard continuous Qi finite elements at a coarse discretization scale H 
as trial functions, whereas the test functions are computed as the solu¬ 
tions of local problems at a finer scale h. The diameter of the support 
of the test functions behaves like mH for some oversampling parameter 
m. Provided m is of the order of log(K) and h is sufficiently small, the 
resulting method is stable and quasi-optimal in the regime where H is 
proportional to In homogeneous (or more general periodic) media, 
the fine scale test functions depend only on local mesh-configurations. 
Therefore, the seemingly high cost for the computation of the test 
functions can be drastically reduced on structured meshes. We present 
numerical experiments in two and three space dimensions. 

Keywords multiscale method, pollution effect, wave propagation, Helmholtz prob¬ 
lem, finite element method 

AMS subject classification 35J05, 65N12, 65N15, 65N30 


1 Introduction 

Standard finite element methods (FEMs) for acoustic wave propagation are 
well known to exhibit the so-called pollution effect IBSOPj . which means 
that the stability and convergence of the scheme require a much smaller 
mesh-size than needed for a meaningful approximation of the wave by finite 
element functions. For an highly oscillatory wave at wave number k, the 
typical requirement for a reasonable representation reads kH < 1 for the 
mesh-size H, that is some fixed number of elements per wave-length. The 
standard Galerkin FEM typically requires at least < 1 where a > 1 

depends on the method and the stability and regularity properties of the 
continuous problem. There have been various attempts to reduce or avoid the 
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Figure 1: Convergence history of the mnltiscale FEM (msPGFEM), the standard Qi 
FEM (FEM) and the best-approximation (bestapprox) in the finite element space 
for a two-dimensional plane wave with wave number k = 2^ (see also Section ^ . 


pollution effect, e.g., discontinuous Galerkin methods [TFObl iFWOfil IFWlll 
IHMP11| . high-order finite elements |MS10l IMSlla] . discontinuous Petrov- 
Galerkin methods |ZMD+1 1 1 [DGMZ12] . or the continuous interior penalty 
method |Wul4) among many others. A good historical overview is provided 


m 


ZMD+11 


The work |Petl4j suggested a multiscale Petrov-Galerkin method for the 
Helmholtz equation where standard finite element trial and test functions 
are modified by a local subscale correction in the spirit of numerical homog¬ 
enization |MP14j . In the numerical experiments of | Pet 14) , a variant of that 
method appeared attractive where only the test functions are modified while 
standard finite element functions are used as trial functions. In this paper, 
we analyze that method and reformulate it as a stabilized Qi method in the 
spirit of the variational multiscale method |Hug95 , |HFMQ9^ IHSOTi IMallli 
IPetlbj . The method employs standard Qi finite element trial functions on 
a grid 9jy with mesh-size H. The test functions are the solutions of local 
problems with respect to a grid 9h at a finer scale h which is chosen fine 
enough to allow for stability of the standard Galerkin FEM over 9fe- The 
diameter of the support of the test functions is proportional to niH for the 
oversampling parameter m. Under the condition that m is logarithmically 
coupled with the wave number n through m ~ log(K), we prove that the 
method is pollution-free, i.e., the resolution condition kH < 1 is sufficient 
for stability and quasi-optimality under fairly general assumptions on the 
stability of the continuous problem. The performance of the method is illus- 
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trated in the convergence history of Figure More detailed descriptions on 
the numerical experiments will be given in Section As the test functions 
only depend on local mesh-configurations, on structured meshes the number 
of test functions to be actually computed is much smaller then the overall 
number of trial and test functions on the coarse scale. In many cases, the 
computational cost is then dominated by the coarse solve and the overhead 
compared with a standard FEM on the same coarse mesh remains propor¬ 
tional to ~ log(K)'^. Even if no structure of the mesh can be exploited to 
reduce the number of patch problems, the method may still be attractive if 
the problem has to be solved many times with different data (same k though) 
in the context of inverse problems or parameter identification problems. 

The paper is structured as follows. Section states the Helmholtz prob¬ 
lem and recalls some important results. The definition of the new Petrov- 
Galerkin method follows in SectionStability and error analysis are carried 
out in Section]^ Section]^ is devoted to numerical experiments. 

Standard notation on complex-valued Lebesgue and Sobolev spaces ap¬ 
plies throughout this paper. The bar indicates complex conjugation and i 
is the imaginary unit. The inner product is denoted by {v,w)i 2 (^q-^ := 
Jq vw dx. The Sobolev space of complex-valued U’ functions over a domain 
u) whose generalized derivatives up to order k belong to L’P is denoted by 
VF^’^(a;;C). The notation A< B abbreviates A < CB for some constant C 
that is independent of the mesh-size, the wave number k, and all further pa¬ 
rameters in the method like the oversampling parameter m or the fine-scale 
mesh-size h] A k, B abbreviates A< B < A. 


2 The Helmholtz Problem 


Let O C for d G {1,2,3}, be an open bounded domain with polyhedral 
Lipschitz boundary which is decomposed into disjoint parts dO. = Tj) U F/j 
with F D closed. The classical Helmholtz equation then reads 


—Au — 

u 

iKU -Vu-v 


f in tl, 

UD on Td, 
9 on Tr 


( 2 . 1 ) 


for the outer unit normal u of H and the real parameter k > 0. For the sake 
of a simple exposition we assume ur = 0. Either of the parts T^ or F/j is 
allowed to be the empty set. In scattering problems, the Dirichlet boundary 
F R) typically refers to the boundary of a bounded sound-soft object whereas 
the Robin boundary T r arises from artificially truncating the full space to 
the bounded domain H |Ihl98) . The variational formulation of (2.1) employs 
the space 

V := W^^(n;C) := {v G : z;|rB = 0}. 
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For any subset a; C we define the norm 

:= + Il^^lli2(a;) any u E y 

and denote ||w||\/ := ||w||y,n- Define on V the following sesquilinear form 

a(v,w) := (Vv,Vw)l2(^) - k^(v,w)i^2(qj - iK(v,w)i^ 2 ^rji)- 

Although the results of this paper hold for a rather general right-hand side 
in the dual of V, we focus on data / E L^(Q;C) and g E L^(r/j;C) for the 
ease of presentation. The weak form of the Helmholtz problem then seeks 
u € V such that 


a{u, v) = if, v)L^{n) + id, v)L2(rn) for ah v € V. 


( 2 . 2 ) 


We assume that the problem is polynomially well-posed |EM12j in the sense 
that there exists some constant 7 (k, D) which depends polynomially on k 
such that 

51?a(r;, w) 


7 (k, D) ^ < inf sup 

^ev\{0}„,ey\{0} llr^llvII'W^II f 


(2.3) 


For instance, in the particular case of pure impedance boundary conditions 
= Fi^, it was proved in |Mel951 ICFOB] by employing a technique of 
|MIB96| that 7 (k, D) < k. Further setups allowing for polynomially well- 
posedness are described in [HetOTl IEM12| IHMP14b] . In particular, the case 
of a medium described by a convex domain (with Robin boundary conditions 
on the outer part of the boundary) and a star-shaped scatterer (with Dirichlet 
boundary conditions) allows for polynomial well-posedness |Het07j . Another 
admissible setting is described in |EM12| where H is a bounded Lipschitz 
domain with pure Robin boundary. For general configurations, however, the 


dependence of the stability constant 7 (k, H) from (2.3) is an open question. 
Throughout this paper we assume that (2.3) is satisfied. The case of a 
possible exponential dependence IBCWG"*"!!) is excluded here. 


3 The Method 

This section introduces the notation on finite element spaces and meshes 
and defines the multiscale Petrov-Galerkin method (msPGFEM) for the 
Helmholtz problem. 

3.1 Meshes and Data Structures 

Let Qh be a regular partition of H into intervals, parallelograms, paral¬ 
lelepipeds for d = 1,2,3, respectively, such that U2h = D and any two dis¬ 
tinct T, T' E Sh are either disjoint or share exactly one lower-dimensional 
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hyper-face (that is a vertex or an edge for d G {2, 3} or a face for d = 3). We 
impose shape-regularity in the sense that the aspect ratio of the elements 
in Sh is uniformly bounded. Since we are considering quadrilaterals (resp. 
hexahedra) with parallel faces, this guarantees the non-degeneracy of the 
elements in Sh- We consider this type of partitions for the sake of a sim¬ 
ple presentation and to exploit the structure to increase the computational 
efficiency. The theory of this paper carries over to simplicial triangulations 
or to more general quadrilateral or hexahedral partitions satisfying suitable 
non-degeneracy conditions or even to meshless methods based on proper 
partitions of unity |HMP14aj . 

Given any subdomain S' C n, define its neighbourhood via 

N(S) :=int(u{rEgH ; rns^0}). 

Furthermore, we introduce for any m >2 the patches 

Ni(S) := N(S) and N”^(S) := N(N"*-i(S)). 

The shape-regularity implies that there is a uniform bound Coi^m = C'oi,m(d) 
on the number of elements in the mth-order patch, 

max cardjiF G Sh '■ K C N™-(T)j < Co\m- 

We abbreviate Co\ ■= Goip. Throughout this paper, we assume that the 
coarse-scale mesh Sh is quasi-uniform. This implies that depends 

polynomially on m. The global mesh-size reads H := max{diam(T) : T G 
9h}- Let Qp{Qh) denote the space of piecewise polynomials of partial degree 
< p. The space of globally continuous piecewise first-order polynomials reads 

The standard Qi finite element space reads 

Vh :=§\SH)nV. 

The set of free vertices (the degrees of freedom) is denoted by 

'■= {z G Q : z is a vertex of and z ^ Td}- 

Let Ih ■ V ^ Vh be a surjective quasi-interpolation operator that acts as a 
stable quasi-local projection in the sense that Ih^Ih = Ih and that for any 
T G Sh and aVi. v GV there holds 

H ^\\v — Ihv\\l^{T) + \\'^^Hv\\l'2(T) ^ (3.1) 

Under the mesh condition that kH < 1 is bounded by a generic constant, 
this implies stability in the || • ||v" norm 

< Cij^y\\v\\v for all v G V, (3.2) 
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with a K-independent constant One possible choice (which we use in 

our implementation of the method) is to define Ih ■= Eh o Hh^ where 11 j;/ is 
the piecewise projection onto Qi{Sh) and Eh is the averaging operator 
that maps Qi{9h) to Vh by assigning to each free vertex the arithmetic 
mean of the corresponding function values of the neighbouring cells, that is, 
for any v £ Qi{9h) and any free vertex z G J^h, 

{Eh{v)){z) = ^ v\t{z) j card{K £^h '■ z£ K}. 

t&Sh ' 

with 2 GT 


Note that EH{v)\ro = 0 by construction. For this choice, the proof of 
(3.11 follows from combining the well-established approximation and stability 
properties of Hh and Eh-, see, e.g., |DE12) . 


3.2 Definition of the Method 

The method is determined by three parameters, namely the coarse-scale 
mesh-size H, and the stabilization parameters h (the fine-scale mesh-size) 
and m (the oversampling parameter) which are explained in the following. 
We assign to any T G Sh its m-th order patch Q-t '■= N™’(T) (for a positive 
integer m) and define for any v,w £ V the localized sesquilinear forms 

anj-iv,w) := - i(Kf, w)L2(r^nanT) 


and 


aT{v,w) := {Vv,Vw)l2^h) - -'i{i^v,w)L2yRndT)- 

Let S/i be a global uniform refinement of the mesh 9 h over H and define 
Vki^r) ■= {v £ Qi{9h) : V = 0 outside Or}. 

Define the null space 


Whi^r) ■= {vh G Vhi^T) ■ Inivk) = 0} 

of the quasi-interpolation operator Ih defined in the previous section. Given 
any nodal basis function G Vh, let £ Wh{V.T) solve the subscale 
corrector problem 


an.j-{w,Xz,T) = aT{w,Az) for all w £ W/j(Dt)- 


(3.3) 


The well-posedness of (3.3) will be proved in Sectionj^ Let Az := X^tsSh 
and define the test function 


A, := A^ - A^. 
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The space of test functions then reads 

Vh '■= span{A^ : z G TT/f}. 

We emphasize that the dimension dimV/f = dimVn is independent of the 
parameters m and h. Figures [2j|^ display typical examples for the test func¬ 
tions and correctors. The multiscale Petrov-Galerkin FEM seeks uh G Vh 
such that 

a{uH,VH) = {f,VH)L^(Q) +{9,VH)L-2(rR) for all UH G E//. (3.4) 

The error analysis and the numerical experiments will show that the 
choice H < k~^, m ^ log(K) suffices to guarantee stability and quasi¬ 
optimality properties, provided that K,°‘h < 1 where a depends on the sta¬ 
bility and regularity of the continuous problem. The conditions on h are the 
same as for the standard Qi FEM on the global fine scale (e.g. < 1 

for stability |Wul4j and < 1 for quasi-optimality |Mel95) in the case of 
pure Robin boundary conditions on a convex domain). 


3.3 Remarks on Generalizations of the Method 

The present approach exploits additional structure in the mesh and thereby 
drastically decreases the cost for the computation of the test functions (A^ : 
z G AIh). Indeed, (3.3) is translation-invariant and, thus, the number of cor¬ 
rector problems to be solved is determined by the number of patch configu¬ 
rations. This number is typically much smaller than the number of elements 
in Shi see Figure]^ for an illustration. 

Some remarks on more general versions of the presented msPGFEM are 
in order. 


Element shapes. As Figurej^illustrates, highly structured meshes are de¬ 
sirable as they lead to a moderate number of patch problems. The method 
presented in Subsection |3.2| considers, for simplicity, a partition of the do¬ 
main in parallelepipeds. While in scattering problems the outer part of the 
boundary T^ results from a truncation of the full space and, hence, the 
choice of a simple geometry (e.g., a cube) is justified, it is extremely im¬ 
portant to guarantee an accurate representation of more general scattering 
objects. This requires more general element shapes such as isoparametric el¬ 
ements or partitions in bricks and simplices with first-order ansatz functions 
on the reference cell (see |Petl4j for simplicial meshes). The msPGFEM and 
its error analysis is also applicable to this situation. The configurations at 
the boundary will then determine the number of corrector problems. 

Fine-scale grid. The present approach is based on a global fine-scale grid 
Sh and a particular choice of the domains V.t, which is convenient for the 
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Figure 2: Coarse-scale trial function (left), corrector Az (middle), and modified 
test function Az = — Xz (right) in ID with k = 2^, H = 2“^, h = 2~^^, m = 2. 



0.8125 0.0625 0.8125 0.0625 


Figure 3: Coarse-scale trial function (left), and element corrector Xz^t (right) in 
2D with K = 2^, H = 2“^, h = 2“^, m = 2 for the patch highlighted in Figured 


implementation of the method. It is, however, not necessary for the domains 
Qt to be aligned with the mesh Also the spaces Wh{^T) can be defined 
over independent fine-scale meshes over VIt- 

Adaptive methods. For certain configurations of the domains VLt) for 
instance in the presence of re-entrant corners, it may be desirable to utilize 
an adaptive fine-scale mesh over for the solution of the corrector problem 
As proven in Lemma below, the corrector problems are coercive 
and mesh-adaptation may improve the efficiency of the fine-scale corrector 
problem. As mentioned in the previous remark, it is indeed possible to 
employ independent fine-scale meshes over different domains Q.j'-, ^k- The 
stability and error analysis for the adaptive case, which are expected to be 
more involved, are not discussed in this paper. 


(3. 
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Figure 4: All possible patch configurations (up to rotations) on a structured mesh 
of the square domain with pure Robin boundary with m = 2. A trial function and 
corresponding corrector for the highlighted patch is depicted in Figure 



Figure 5: Modulus of the idealized test function for m = oo, H = 2“^, h = 2“^ 
in 2D in a logarithmically scaled plot. The dots indicate the grid points of the 
coarse mesh. Left: k = 2®; right: k = 2®. 
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4 Error Analysis 

We denote the global finite element space on the fine scale by 14 := 14(0) = 
{9h)GV. We denote the solution operator of the element corrector problem 
( |3.3| ) by QT,m- Then any z G 'Nh and any T G 9h satisfy Xz,t = QT,m{^z) 
and we refer to Ct,™ as element correction operator. The map e->■ Az 

described in Subsection 3.2 defines a linear operator Qm via Cm(A^) = 
for any z G ISl//, referred to as correction operator. For the analysis we in¬ 
troduce idealized counterparts of these correction operators where the patch 
Q.T equals 12. Define the null space space M4 '■= {v ^Vh '■ Ih{v) = 0}. For 
any t; G lA, the idealized element corrector problem seeks Qt,ooV G 114 such 
that 

a{w, Gt,oo'v) = ariwjv) for all w G 114. (4.1) 

Furthermore, define 

GooV := ^T,ooV. (4.2) 

t&9h 

It is proved in [MSlOl Corollary 3.2] that the form a is continuous and 
there is a constant Ca such that 

a{v,w) < Calif ||y||r(;||\/ for all v,w £V. 

The following result implies the well-posedness of the idealized corrector 
problems. 

Lemma 1 (well-posedness of the idealized corrector problems). Provided 

Ci„y^lHK<l/V2, (4.3) 

we have for all w G 114 equivalence of norms 

l|Vw;||i2(n) < |h||y < vW2||Vf;|42(Q) 

and ellipticity 

< ^a{w,w). 

Proof. For any w G 114 the property ( |3.1[ ) implies 

k2||(1 - Ih)w\\1,^^^ < □ 

Lemma implies that the idealized corrector problems (4.2) are well- 
posed and the correction operator Soo is continuous in the sense that 

||Coo'y//||y < Ce||f_H-||y for all vh £ 14 

for some constant Cq ~ 1. Since the inclusion Whi^lr) 4 114 holds, the well- 
posedness result of Lemma carries over to the corrector problems (3.3) in 
the subspace Wh{^T) with the sesquilinear form a^j.. 
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The proof of well-posedness of the Petrov-Galerkin method (3.4) will be 


based on the fact that the difference (Coo ~ Qm){v) decays exponentially with 
the distance from supp(u). In the next theorem, we quantify the difference 
between the idealized and the discrete correctors. The proof will be given 
in Appendix of this paper and is based on the exponential decay of the 
corrector C 00 A 2 itself, see Figure [5 That figure also illustrates that the decay 
requires the resolution condition (4.3), namely kH < 1. 


Theorem 1. Under the resolution condition (4.3) there exist constants Ci ~ 
1 ~ C 2 and 0 < < 1 such that any v G Vh, any T G 9h and any m G N 

satisfy 


||V(eT,oo^^ - eT,mv)\\mn) < Ci/3-||Vu||i2(r), (4.4) 

||V(eoo^; - emv)\\mn) < C2^/C^^\\^v\\L^n)■ (4.5) 

□ 


Provided h is chosen fine enough, the standard FEM over S/x is stable 
in the sense that there exists a constant Cfem such that with 7 (k, n) from 
(2.3) there holds 


(C'fem 7(«^) ^ < inf sup 


(Ra(u, w) 




(4.6) 


This is actually a condition on the fine-scale parameter h. In general, the 
requirements on h depend on the stability of the continuous problem |Mel95j . 

Theorem 2 (well-posedness of the discrete problem). Under the resolution 
conditions ( |4.3[ ) and (4.6) and the following oversampling condition 

m> |log {V^Cayf^lCijjCIjj yC 2 ^ Co\,mC-F'EM'y{l^, f^)) | / |log(/3)|, (4.7) 


problem (3.4) is well-posed and the constant Cpq := 2C'7^yC'eC'FEM satisfies 

iRa{vH,VH) 


(C'pg 7(^) ^)) ^ < inf sup 


\\vh\\v\\vh\\v 


Proof. Let uh G Vh with IluFflly = 1. From (4.6) we infer that there exists 
some V £ Vh with ||u||f = 1 such that 


ika{uH — QooiuH),v) > (C'fEM7(^) ^)) ^W'^H — QooiuH)\\v- 


It follows from the structure of the sesquilinear form a that Coaiun) solves 
the following adjoint corrector problem 

a{Qcxi{'dH)iW) = a{uH,w) for all re G VF/i, (4.8) 
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cf. |MSllbl Lemma 3.1], Let vh '■= (1 — Qm)lHV G Vh- We have 

a{uH, vh) = a{uH, (1 - Qoo)Ihv) + a{uH, (Coo - Qm)lHv). (4.9) 

Since Coo is a projection onto VF/j, we have (1 — Coo)(l — Ih)v = 0 and, thus, 
(1 — Qoo)Ihv = (1 — Coo)?^- The solution properties (4.8) of Qoo{uh) and 
([4d|-([4(^ of QooV prove a{uH,Coov) = a{QooiuH),v). Hence, 


?R.a{uH, (1 - Qoo)Ihv) = ?R.a{uH - Coo{uh), v) 

> (C'fem7(«^)^)) ^\\uh - Qoo{uh)\\v■ 
Furthermore, the estimate (|3.2|) implies 


1 = Ik^llv = WlniuH - Coo('Uh))||v < Cij^y\\uH - Qod{uh)\\v- 


The second term on the right-hand side of (4.9) satisfies with ||riH||v' = 1 
and Lemma [T] that 

\a{uH, (Coo — Qm)lHv)\ < -\/3/2Ca||V(Coo “ ||L2(O)• 

Altogether, it follows that 


‘?R.a{uH,VH) > 


y V(Coo ^m)lHv\\L2(^Q'j 


CijjyCpEM'yiK, H) 

Theorem]^ and (3.1) show that 

||V(Coo - en.)lHv\\L^(^n) < C2,/a^Pn\VlHv\\ < v^/ 3 ™. 

— Coo)'y||y < Ce imply the 

□ 


Hence, the condition ( |4.7[ ) and H'ChIIv = 
assertion. 


Remark 1 (adjoint problem). Under the assumptions of Theorem^ prob¬ 
lem (3.4) is well-posed and, thus, it follows from a dimension argument that 
there is non-degeneracy of the sesquilinear form a over Vh x Vh- Thus, the 
adjoint problem to (3.4) is well-posed with the same stability constant as in 
Theorem 

The quasi-optimality result requires the following additional condition 
on the oversampling parameter m. 


m > |log { 202 ^Tl 


|log(/3)|. 


(4.10) 


Theorem 3 (quasi-optimality). The resolution conditions (4.3) and (4.6) 
and the oversampling conditions (4.7) and (4.10) imply that the solution uh 
to (3.4) with parameters H, h, and m and the solution Uh of the standard 
Galerkin FEM on the mesh Sfe satisfy 


\uh - uhWv - lH)uh\\v ^ min \\uh-VH\\v- 

vh&Vh 
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Proof. Let e := Uh — ujj. The triangle inequality and Lemma yield 
||e||v/ < 11(1 - lH)uh\\v + \\lHe\\v- 


It remains to bound the second term on the right-hand side. The proof 
employs a standard duality argument, the stability of the idealized method 
and the fact that our practical method is a perturbation of that ideal method. 
Let zh £ Vh be the solution to the dual problem 

V/^e) -I- = a{vH, (1 - Qoo)zh) 

for all Vh £ Vh (cf. Remark Q. The choice of the test function vh = Ih^ 
implies that 


(1 Gqc)zh) — (6^ Coo).2/f)(1 Qm)zH\ 


The identity lHiQm — Qoo)zH = 0, the resolution condition (4.3), the estimate 
( |4.5[ ), and the stability of the adjoint problem imply for the first term on the 
right-hand side that 


(Cm ^oo^Zh') 

< Ca^/sJ2\\IHe\\v\\V{Qm - Soo)^h||l2(q) 

< C2^/C^Cay/W^\\lHe\\vr\\^ZH\\ 

< C2,/C^C^aCpGli^,n)y/W^r\\lHefv 


The condition ( 4.10| ) implies that this is < 2 l|L//e||y. The Galerkin orthog¬ 
onality a{uh — Uh, (1 — Qm)ZH) = 0, the solution property (4.2) of QooZh, 
the resolution condition ( |4.3| ) and the exponential decay (4.5) imply for the 
second term 


u{Ih(^, (1 ^m)ZH^ — u(^IhU}i Uh, (1 Gm)ZH) 

— u{lHUh Uh, (Coo ^m)ZH') 

< \/3/2C'aC'2y^C'ol,m/3™'||L_ffn/i — Uh\\v\\'VZ hWl^^Q)- 
The stability of the adjoint problem implies 


Thus, 


\\'^ZH\\L2{n) < CpG'y{K,^)Ca\\lHe\\v- 


a{lHe, {l — Qm)zH) < \/3/2CQC'2-y/C'ol,mC'PG/3™'7(^) ^)\\lHUh—Uh\\v\\lHe\\v- 


The term ||///e||y can be absorbed and the oversampling condition (4.7) 
implies that ■\/Coi^'y{K, Q) is controlled by some ^-independent constant. 
The combination with the foregoing displayed formulae concludes the proof. 

□ 


13 



















Multiscale Petrov-Galerkin FEM 


The following consequence of Theorem [^states an estimate for the error 
u — uh- 

Corollary 1. Under the conditions of Theorem^ the discrete solution uh 
to (3.4) satisfies with some constant C ~ 1 that 


\u — uh\\v ^ \\u — Uh\\v + C min \\uh — VH\\v- 

vh&Vh 


In particular, provided that the solution satisfies u G for 0 < s < 1, 

the error decays as ||n — □ 

Remark 2. In the idealized case that m = oo, we have Uh — InUh £ Wh 
and, thus. 


a{uh - Ihuk, (1 - Qoo)vh) = 0 for all vh G Vh- 
Therefore, problem (|3.4[) and the Galerkin property show that uh = IhUh- 


5 Numerical Experiments 

We investigate the method in three numerical experiments. The convergence 
history plots display the absolute error in the norm || • ||y versus the mesh 
size H. 


5.1 Plane Wave on the Square Domain 

On the unit square Q = (0,1)^, we consider the pure Robin problem Tr = dQ 
with data given by the plane wave u(x) = exp (—zkx • ( qJ )). 

Figure 6a-6c displays the convergence history for k = 2®, 2'^, 2^ and the 
fine-scale mesh parameter h = 2“^^. The best-approximation error of con¬ 
tinuous Qi functions in || • ||y and the error of the standard Galerkin FEM on 
the same coarse mesh are plotted for comparison. As expected, the standard 
FEM clearly exhibits the pollution effect, and larger values of k increase the 
discrepancy between the approximation error of the FEM and the theoret¬ 
ical best-approximation by Qi functions in the regime under consideration. 
In contrast, the approximation by the msPGFEM can compete with the 
best-approximation on meshes that allow a meaningful representation of the 
solution. We stress the fact that the convergence history plots merely take 
into account the coarse mesh-size H, but the computational cost in the mul¬ 
tiscale method is moderately higher than in the standard FEM due to the 
increased communication caused by the coupling m ~ log(K). 

For the oversampling parameter m = 2, the number of corrector problems 
to be solved for the finest mesh 9 h is 49 out of 1 048 576 when no symmetry 
is exploited. 
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Figure 6: (a)-(c). Comparison of the msPGFEM with the best approximation in 
11 • IIV and the standard Galerkin FEM for the 2D plane wave example for n = 
2®, 2^, 2®. (d). Dependence on the fine mesh parameter h in the 2D plane wave 
example with k = 2®. 
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Figure 7: Convergence history for the 2D plane wave example with k = 2^ (left) 
and K = 2®(right), h = and varying m. 


Figure 6d displays the dependence on the fine mesh parameter h for 
K = 2® and oversampling parameter m = 6. Since the multiscale method 
based on the fine grid Sh computes approximations of the FEM solution on 
that fine grid, e.g. uh = InUh for m = oo as in Remark it is clear that the 
accuracy of the msPGFEM is limited by the accuracy of the standard FEM 
on the fine scale. This can be observed in Figure [6d| It can be also seen that a 
finer fine-scale mesh-size h improves the error of the msPGFEM towards the 
best-approximation. In this two-dimensional example, the quasi-optimality 
constant appears to be close to 1 

Next, we study the dependence on the oversampling parameter m. Fig¬ 
ure displays the convergence history for k, = 2^ and k = 2®. The fine mesh 
parameter is /i = 2“^^ and m varies from m = 1 to m = 6. It turns out 
that for the present configuration, the value m = 2 is sufficient for quasi¬ 
optimality. In the range where H is significantly larger than k~^ and the 
resolution condition is violated, larger oversampling parameters may lead to 
larger errors, which is not surprising in view of the lack of decay, see also 
Figure]^ This, however, is no more the case as soon as H is small enough 
to allow for a meaningful representation of the wave. 


5.2 Multiple Sound-Soft Scatterers in 2D 

We consider the domain 


^ - (0> 1)^( X [^, ^] U [i, if] X [A, u [A, A] X [A J] ) 

from Figure]^ The incident wave u\yi{x) = exp(—i/tx • (o'g)) is incorporated 
through the Robin boundary condition with g := ZKUin -|- di/Uin on the outer 
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Figure 8: 


tion 5.2 


Coarse mesh for the square domain with three scatterers from Subsec¬ 


boundary Fjj := {x G {0,1} or y G {0,1}}. On the remaining part of 
the boundary F/r := dQ \ F/j we impose homogeneous Dirichlet conditions. 
We choose the hne mesh parameter as h = 2~^^. Since the exact solution 
is unknown, we compute a reference solution with the standard Qi FEM 
on the hne mesh S/j and we compare the coarse approximation with this 
reference solution. Errors committed by the hne scale are not included in the 
discussion. Figure [^displays the convergence history for k = 2® and k = 2®. 
The oversampling parameter m varies from m = 1 to m = 4. As in the 
foregoing example, the value m = 2 for the oversampling parameter seems 
to be sufficient for the quasi-optimality and even a quasi-optimality constant 
close to 1 in the range of wave numbers considered here. In particular, the 
pollution effect that is visible for the standard Galerkin FEM is not present 
for the msPGFEM. Reduced convergences rates which are expected from the 
presence of re-entrant corners are not visible in this computational range. For 
the oversampling parameter m = 2, the number of corrector problems to be 
solved for the hnest mesh 9 h is 210 out of 61952 when no symmetry is 
exploited. 


5.3 Plane Wave on the Cube Domain 


On the unit cube P = (0,1)^, we consider the pure Robin problem with data 
given by the plane wave u{x) = exp(-zKx • ^ ). 

We choose k = 2^. Figure 10 compares the error of the msPGGEM 
h = 2~^ and m G {1,2,3,4} with the best-approximation in the || • ||y 
norm and the error of the standard Galerkin FEM. Also in this example, 
the msPGFEM is pollution-free for the oversampling parameter m > 2. 
The quasi-optimality constant appears slightly larger than in 2D. For the 
oversampling parameter m = 2, the number of corrector problems to be 
solved for the hnest mesh 9 h is 343 out of 262 144 when no symmetry is 
exploited. 
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Figure 9: Convergence history for the multiple scattering example from Subsec- 
for K = 2^ (left) and k = 2® (right) h = 2~^^. 


tion 


5.2 



Figure 10: Convergence history for the 3D plane wave example for k = 2® and 
h = 2-\ 
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A Appendix: Proof of Theorem 

For the sake of completeness we also present a proof of the exponential decay 
result Theorem which is central for the method. The idea of the proof is 
the same as in the previous proofs of the exponential decay |MP141 IHP131 
lEGMPl^ IHMP14al IBP 14] in the context of diffusion problems. The differ¬ 
ence especially with respect to |Petl4j is that here the quasi-interpolation is 
a projection. This simplihes the proofs and leads to slightly better rates in 
the exponential decay that have been experimentally observed in |Petl4) . 

Let J/j : (7^(0) —)• Vh denote the nodal Qi interpolation operator. Stan¬ 
dard interpolation estimates and the inverse inequality prove for any T G Sh 
and all q G Q 2 {T) the stability estimate 

\\^Ihq\\L‘^(T) < C'/hll'^^llL2(T)- A-l) 

In the proofs we will frequently make use of cut-off functions. We collect 
some properties in the following lemma. 

Lemma 2. Let rj G §^{5h) be a funetion with values in the interval [0,1] 
satisfying the bound 

W^vh^iu) < Cr,H-^ (A.2) 

and let Ik := supp(Vr 7 ). Given any subset JC C Sh, any cj) G satisfies for 
S = UJC C n that 


UWl^s) < ^||V0|U2(n(s)) (A.3) 

11(1 - lH)Ih{vmLfiS) < (A.4) 

\\'^iv4‘)\\L^{S) ^ ll^<('llL2(5n{supp(»?)}) + (A.5) 

Proof. The property ( |3.1| ) readily implies ( |A.3| ). Furthermore, ( |3.1| ) implies 

11(1 - lH)Ih{m\LfiS) < HCj,V^y\\S/hivmLfiN(S)y 


Estimate ( A.lj) leads to 


l|V//i(r/(())||L2(M(5)) < C'4||V(?7(?i)||i2(M(5)). 

This proves ( |A.4| ). For the proof of (A.5) the product rule and (A.2) imply 

II^(^</’)IIl 2(5) < ll^<(’llL2(Sn{supp(77)}) + ^ll</’llL2(5nS)- 

The combination with (|A.3|) concludes the proof. □ 


Theorem 4 (decay). Under the resolution eondition (4.3), there exists 0 < 
fi < 1 sueh that, for any vh £ Vh cLnd all T G Sh and m G N, 


||VeT,oo'yH||L2(Q\N"*(T)) < CI3'^\\\7vh\\l^{T)- 
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Proof. We define the cut-off function rj £ §^{Sh) via 

r] = 0 and rj = I inO\N™-2(r). 

Note that r] is thereby also uniquely defined on the set 31 := supp(Vr/). The 
shape-regularity implies that rj satisfies (A.2). Let vh £ Vh and denote 
(f := Qt,ooVh £ Wh- Elementary estimates lead to 

ll^'^lln\N-(r) ^ l(V</>,r/V(/))i2(Q)| < |(V(/>, V(r?(/>))i2(f^)| -F |(V(/>, (/>V77)i2(Q)| 

< Mi + M 2 + M 3 + M 4 ^ 


for 

Ml :=|(V</),V((l-4)(r?</))))i2(o)| M 2 := - lH)Ihiv4>)))LHn)\ 

M 3 := \{V(l),V{lHlh{v(l>)))L^(n)\ M 4 := |(V(/>, (/>V? 7 )i 2 (Q)|. 

The property ( |A.l ) proves 

Ml < ||V(/>||i2(3j) ||V(#-4(77(/.))||i2(3j) < ||V(/>||i;,2(3j)||V(#)||i2(3j). 


Hence, it follows with (A.5) that 


^1 ^ I|v</’IIl2(d?)I|v</’IIl2(n(3j))- 


Since w := {1 — lH)Ihiv4') ^ the identity (4.1) and the fact that the 
support of w lies outside T imply a{w,4>) = aT{w,VH) = 0 and therefore 

M 2 = a{w,(l)) -F k ^{ w , 4 >) = K^iw.cj)) < K^||w||L2(N(S))ll0llL2(N(a?))- 


The estimates (A.3) and (A.4) and the resolution condition kH < 1 from 
( 4.3 ) imply 

M 2 < ||V(/>||x,2(M2(3j))||V(r?0)||i;^2(M2(3j)). 

The application of ( |A.5 ) yields 

M 2 < ||V(/>||2,2(N2(3J))(||V(/>||i^2(M2(3J)) -F ||V0||2,2(N(3J))) < ||V(?i||^2(M2(3^)). 

The function lHlh{'>l4’) vanishes outside N(3J). Hence, the stability and ap¬ 
proximation properties ( |3.1[ ) and ( |A.l ) lead to 

M 3 < ||V(/)||2,2(n(3j))||V(/hI/i(770))||i;^2(m(3j)) 

^ ll^'/’llL2(N(3?))||V(r/(/))||2,2(N2(3j)). 

With ( A.5| ) we obtain 

M 3 < ||V(/>||i;^2(m(3j))(||V(/)||2,2(n 2(3^)) -I- ||V(/>||i;,2(N(3i))) < ||V0||^2 (n2(3j))• 
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For the term M 4 , the Lipschitz bound ( |A.2[ ) and (A.3) prove 

Altogether, it follows for some constant C that 

ll'^‘^llL2(f2\N’"(T)) - ^II^'?^IIl2(n2(s))- 
Recall that N2(3J) = N™(r) \ N"^-5(r). Since 

ll^'^llL2(r2\N™(r)) + ll^'/’lli2(N ’"(T)\N"*-5(T)) ~ ll^</’llL 2 (n\N— 5(r))> 
we obtain 


+ c )||V(?i||^2(f7\M™(r)) — ll^‘?^llL2(f7\N'"-5(T))- 


The repeated application of this argument proves for /3:=(1 + C' ^<1 
that 

This is the assertion. □ 

We proceed with the proof of Theorem 

Proof of Theorem^ We define the cut-off function rj G §>^{9h) via 

r] = 0 in \ N"*-^(r) and r/= 1 in N™- 2 (r). 

This function is thereby uniquely defined and satisfies the bound ( |A.2 ). Since 
(1 — lH)Ih{'n^T,oov) G Whi^T)-, we deduce with Cea’s Lemma, the identity 
IhQt, coV = 0 and the approximation and stability properties (3.1) and (A.l) 
and the resolution condition (4.3) that 

||V(er,oo'y - CT,m^')|li2(^7) < ||Sr,oo'y - (1 - lH)IhivQT,ocV)\\v 

= 11(1 - lH)Ih{QT,ooV - ??Cr,oot’)||y,n\{r,=l} 

< ||V(1 - 

Note that N(n \ {r/ = 1}) = \ N™'“^(T). This and Theorem 0 prove ( |4.4[ ). 

Define z := (Coo — Qm)v and zt '■= {Qt,oo — QT,m)v. The ellipticity from 
Lemma proves 


JI^^IIl 2 (q) < 


^ a{z,ZT) 

T&Sh 
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We define the cut-off function rj G S^(Sh) via 

r]=l mn\N'^+^{T) and rj = 0 in N™+^(r). 

This function is thereby uniquely defined and satisfies the bound ( |A.2 ). For 
any T G we have (1 — lH)Ih{fl^) ^ with support outside VIt- Hence, 
we obtain with z = I^z that 

a(z, zt) = a{Ih{z - rjz), zt) + a^Hlhivz), zt). 

The function z — vanishes on S '.= {r] = 1}. Hence, the first term on 

the right-hand side satisfies 

\a{Ih{z - r]z),ZT)\ < Ca\\Ih{z - Vz)\\v,n\s\\zT\\v- 

The Friedrichs inequality with constant Cp proves together with the stability 
( A.l I and the estimate ( |A.5 ) applied to the cut-off function (1 — rj) that 

Whiz - 'nz)Wv,n\s ^ Vl + (CfkH')2||Vz||2,2(q\ 5) < ||Vz||x,2(pj\5). 

Furthermore, Inhivh vanishes on H \ N(supp(l — r])). Hence, we infer 
from Friedrichs’ inequality and the resolution condition (4.3), the stability 
properties (3.1) and (A.l) and the ( |A.5[ ) that 

\(^izT 1 1 Hhivh}\ ~ II IIl^(N2(supp(1—77))) IIII y• 

The sum over all T G 9h and the Cauchy inequality yield with the finite 
overlap of patches 

iiv^iii.,„, < E II II (N2(supp(l— 77))) II II y 
TGSh 


\/C'ol,777|| V^||_L^(n) . / II'^'TI 

V TGSh 


2 

V 


The combination with (4.4) concludes the proof. 


□ 


References 

[BCWG+ll] T. Betcke, S. N. Chandler-Wilde, I. G. Graham, S. Langdon, and M. Lind¬ 
ner. Condition number estimates for combined potential integral operators 
in acoustics and their boundary element discretisation. Numer. Methods 
Partial Differential Equations, 27(l):31-69, 2011. 

[BP14] D. Brown and D. Peterseim. A multiscale method for porous microstruc¬ 

tures. ArXiv e-prints, November 2014. 

[BSOO] I. M. Babuska and S. A. Sauter. Is the pollution effect of the FEM avoidable 

for the Helmholtz equation considering high wave numbers? SIAM Rev., 
42(3):451-484, 2000. 


22 


















Multiscale Petrov-Galerkin FEM 


[CF06] 

[DE12] 

[DGMZ12] 

[EGMP13] 

[EM12] 

[FW09] 

[FWll] 

[Het07] 

[HFMQ98] 

[HMPll] 

[HMP14a] 

[HMP14b] 

[HP 13] 
[HS07] 

[Hug95] 

[Ihl98| 

[Main] 


P. Gummings and X. Feng. Sharp regularity coefficient estimates for 
complex-valued acoustic and elastic Helmholtz equations. Math. Models 
Methods Appl. Set., 16(1):139-160, 2006. 

D. A. Di Pietro and A. Ern. Mathematical aspeets of discontinuous Galerkin 
methods, volume 69 of Mathematiques & Applications (Berlin). Springer, 
Heidelberg, 2012. 

L. Demkowicz, J. Gopalakrishnan, I. Muga, and J. Zitelli. Wavenumber ex¬ 
plicit analysis of a DPG method for the multidimensional Helmholtz equa¬ 
tion. Comput. Methods Appl. Mech. Engrg., 213/216:126-138, 2012. 

D. Elfverson, E. H. Georgoulis, A. Malqvist, and D. Peterseim. Convergence 
of a discontinuous Galerkin multiscale method. SIAM Journal on Numerical 
Analysis, 51(6):3351-3372, 2013. 

S. Esterhazy and J. M. Melenk. On stability of discretizations of the 
Helmholtz equation. In Numerical analysis of multiscale problems, volume 83 
of Leet. Notes Comput. Sei. Eng., pages 285-324. Springer, Heidelberg, 2012. 

X. Feng and H. Wu. Discontinuous Galerkin methods for the Helmholtz 
equation with large wave number. SIAM J. Numer. Anal, 47(4):2872-2896, 
2009. 

X. Feng and H. Wu. /ip-discontinuous Galerkin methods for the Helmholtz 
equation with large wave number. Math. Comp., 80(276): 1997-2024, 2011. 

U. Hetmaniuk. Stability estimates for a class of Helmholtz problems. Com- 
mun. Math. Set., 5(3):665-678, 2007. 

T. J. R. Hughes, G. R. Feijoo, L. Mazzei, and J.-B. Quincy. The varia¬ 
tional multiscale method—a paradigm for computational mechanics. Com¬ 
put. Methods Appl. Mech. Engrg., 166(l-2):3-24, 1998. 

R. Hiptmair, A. Moiola, and I. Perugia. Plane wave discontinuous Galerkin 
methods for the 2D Helmholtz equation: analysis of the p-version. SIAM J. 
Numer. Anal, 49(l):264-284, 2011. 

P. Henning, P. Morgenstern, and D. Peterseim. Multiscale Partition of Unity. 
In M. Griebel and M. A. Schweitzer, editors, Meshfree Methods for Partial 
Differential Equations VII, volume 100 of Leeture Notes in Computational 
Science and Engineering. Springer, 2014. 

R. Hiptmair, A. Moiola, and I. Perugia. Trefftz discontinuous Galerkin 
methods for acoustic scattering on locally refined meshes. Appl. Numer. 
Math., 79:79-91, 2014. 

P. Henning and D. Peterseim. Oversampling for the multiscale finite element 
method. Multiscale Modeling & Simulation, 11(4):1149-1175, 2013. 

T. Hughes and G. Sangalli. Variational multiscale analysis: the fine-scale 
Green’s function, projection, optimization, localization, and stabilized meth¬ 
ods. SIAM J. Numer. Anal., 45(2):539-557, 2007. 

T. J. R. Hughes. Multiscale phenomena: Green’s functions, the Dirichlet- 
to-Neumann formulation, subgrid scale models, bubbles and the origins of 
stabilized methods. Comput. Methods Appl. Mech. Engrg., 127(l-4):387-401, 
1995. 

Frank Ihlenburg. Finite element analysis of acoustic scattering, volume 132 
of Applied Mathematical Sciences. Springer-Verlag, New York, 1998. 

A. Malqvist. Multiscale methods for elliptic problems. Multiscale Model. 
SimuL, 9(3): 1064-1086, 2011. 


23 



Multiscale Petrov-Galerkin FEM 


[Mel95] 

[MIB96] 

[MP14] 

[MSIO] 

[MSlla] 

[MSllb] 

[Pet 14] 
[Pet 15] 
[TF06] 

[Wul4] 

[ZMD+11] 


J. M. Melenk. On generalized finite-element methods. ProQuest LLC, Ann 
Arbor, MI, 1995. Thesis (Ph.D.)-University of Maryland, College Park. 

Ch. Makridakis, F. Ihlenburg, and I. Babuska. Analysis and finite element 
methods for a fluid-solid interaction problem in one dimension. Mathematical 
Models and Methods in Applied Sciences, 06(08):1119-1141, 1996. 

A. Malqvist and D. Peterseim. Localization of elliptic multiscale problems. 
Math. Comp., 83(290):2583-2603, 2014. 

J. M. Melenk and S. A. Sauter. Convergence analysis for finite element dis¬ 
cretizations of the Helmholtz equation with Dirichlet-to-Neumann boundary 
conditions. Math. Comp., 79(272):1871-1914, 2010. 

J. M. Melenk and S. Sauter. Wavenumber explicit convergence analysis for 
Galerkin discretizations of the Helmholtz equation. SIAM J. Numer. Anal., 
49(3):1210-1243, 2011. 

J. M. Melenk and S. A. Sauter. Wave-number explicit convergence analysis 
for Galerkin discretizations of the Helmholtz equation. SIAM J. Numer. 
Anal., 49:1210-1243, 2011. 

D. Peterseim. Eliminating the pollution effect in Helmholtz problems by 
local subscale correction. ArXiv e-prints, 1411.1944, 2014. 

D. Peterseim. Variational multiscale stabilization and the exponential decay 
of fine-scale correctors. ArXiv e-prints, 1505.07611, 2015. 

R. Tezaur and C. Farhat. Three-dimensional discontinuous Galerkin el¬ 
ements with plane waves and Lagrange multipliers for the solution of 
mid-frequency Helmholtz problems. Intemat. J. Numer. Methods Engrg., 
66(5):796-815, 2006. 

H. Wu. Pre-asymptotic error analysis of CIP-FEM and FEM for the 
Helmholtz equation with high wave number. Part I: linear version. IMA 
J. Numer. Anal, 34(3): 1266-1288, 2014. 

J. Zitelli, I. Muga, L. Demkowicz, J. Gopalakrishnan, D. Pardo, and V.M. 
Calo. A class of discontinuous Petrov-Galerkin methods, part IV: The op¬ 
timal test norm and time-harmonic wave propagation in ID. Journal of 
Computational Physics, 230(7):2406 - 2432, 2011. 


24 



